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The 60th birthday of Johann Rafclski was celebrated during the Strange 
-ness in Quark Matter 2011 in Krakow. Johann was born in Krakow and 
he initiated the series of the SQM conferences. This report, which briefly 
presents my personal view on a history of multi-particle production in high 
energy collisions, is dedicated to Johann. 

PACS numbers: 25.75.-q, 25.75.Ag, 25.75.Nq 

1. Introduction 

The systematic study of particle production in collisions of high energy 
particles started about 60 years ago with the construction of first accelera- 
tors. Here, I will briefly present my personal view on the history of this era. 
Johann Rafelski, whose 60th birthday was celebrated at the Strangeness in 
Quark Matter conference in Krakow has been one of the key contributors 
since the mid- 70s [U [2]. Therefore his impact on the development of the 
field will be emphasized. 



* Presented at Strangeness in Quark Matter 2011, Krakow, Poland 
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2. Experimental and theoretical status quo 



The experimental and theoretical status quo of multi-particle produc- 
tion in high energy collisions is summarized in Fig. [TJ where a sketch of the 
transverse mass, my, spectra of hadrons produced in p+p interactions at 
the center-of-mass energy y/s = 50 GeV is shown [3]. Clearly, there are 
three distinct domains: 

1. the soft domain, < 2 GeV, in which a predominant majority of 
all particles is produced, tut spectra are exponential and produced 
particles are essentially uncorrelated; their production properties are 
well described by statistical and hydrodynamical models, 

2. the hard domain, 2 < my < \fs/2 GeV, in which only a small fraction 
of all produced particles is located, spectra follow a power-law 
dependence and produced particles are strongly correlated (organized 
in so-called jets); in this domain particle production properties are 
best described by dynamical QCD-based approaches, 

3. the threshold domain, tut ~ \fs/2, the spectrum in this domain is 
not measured due to a very low particle yield, it is believed to steeply 
decrease to zero with mx approaching its threshold value given by 
energy-momentum conservation laws, \fs/2. 




SKETCH: p+p at 50 GeV 
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Fig. 1. Experimental and theoretical status quo of multi-particle production in 
high energy collisions. 
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3. History of mult i particle production in short 

A brief history of multi-particle production is shown in Fig. [2j The field 
started in the 50s with discoveries of hadrons, first in cosmic-ray experi- 
ments and soon after in experiments using beams of particles produced in 
accelerators. Two classes of models were developed in this time, namely sta- 
tistical and dynamical models of hadron production. The latter class was 
initiated by the scattering matrix approach. Systematic results on prop- 
erties of hadrons and their interactions with electrons accumulated in the 
60- 70s led to discoveries of sub-hadronic particles, quarks and gluons. The 
idea of the quark-gluon plasma (QGP) was formulated. Subsequently, sta- 
tistical approaches of quark and gluon hadronization and a statistical model 
of quark and gluon production were developed. In parallel, a dynamical the- 
ory of strong interactions (QCD) was established. Many QCD-based and 
QCD-inspired models of multi-particle production appeared. Finally, ex- 
perimental studies of nucleus-nucleus collisions at high energies resulted in 
discoveries of strongly interacting matter and its phase transition between 
its hadron gas and quark-gluon plasma forms. 

Johann Rafelski greatly contributed to the field by developing both sta- 
tistical and QCD-inspired approaches to the quark-gluon plasma and its 
hadronization. For a long time his ideas inspired experimental developments 
and he was always devoted to the interpretation of new data. 
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Fig. 2. History of multi-particle production in short. Johann Rafclski's contribu- 
tions span between statistical and QCD-inspired approaches and are deeply rooted 
in experiment. 
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4. Discoveries of hadrons 

Naturally, the first hadrons, discovered in collisions of cosmic-ray par- 
ticles, were the lightest ones, pion, kaon and A. With the rapid advent of 
particle accelerators (see Fig. [3] for a brief history) new particles were un- 
covered almost day-by-day. There are about 1000 hadronic states known 
so far. Their density in mass increase approximately exponentially [1] as 
predicted by Hagedorn's Statistical Bootstrap Model [5]. 

Pioneering discoveries with cosmic-rays: 

-1947: pion (emulsion, Powell et al.) 

-1947: kaon and A (cloud chamber, Rochester, Butler) 



Systematic studies with accelerators, 
in particular: 



-1953: 


Cosmotron at BNL - 


3 GeV 


-1954: 


Bevatron at LBL - 


3 GeV 


-1959: 


PS at CERN - 


28 GeV 


-1960: 


AGS at BNL 


33 GeV 


-1967 


U-70 at IHEP 


70 GeV 


-1976: 


Main Ring at FNAL 


500 GeV 


-1976: 


SPS at CERN 


400 GeV 


-2000: 


RHIC 20 000 GeV 


-2009: 


LHC 30 000 000 GeV 



2010: about 1000 
hadronic states 




Fig. 3. Discoveries of hadrons and a brief history of particle accelerators. The max- 
imum beam energy of accelerators is given in the fixed target system. Accelerators 
used so far in the study of multi-particle production are indicated in red. 
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5. Statistical hadron production 

The first statistical model of multi-hadron production was proposed by 
Fermi [6]. He assumed that hadrons produced in high energy collisions 
are in equilibrium and that the energy density of the created hadronic sys- 
tem increases with increasing collision energy. Soon after, Pomeranchuk [7] 
pointed out that hadrons cannot decouple (freeze-out) at high energy densi- 
ties. They will rather continue to interact while expanding until the matter 
density is low enough for interactions to be neglected. He estimated the 
freeze-out temperature to be close to pion mass, ~150 MeV. Inspired by 
this idea Landau [8], and his collaborators formulated a quantitative hydro- 
dynamical model describing the expansion of strongly interacting hadronic 
matter between Fermi's equilibrium high density stage (the early stage) and 
Pomeranchuk's low density decoupling stage (the freeze-out). The Fermi- 
Pomeranchuk-Landau picture serves as a base for modeling high energy 
nuclear collisions up to now [S]. 

In the 60s Hagedorn made an important conjecture, namely that matter 
composed of hadrons has a maximum temperature, the so-called Hagedorn 
temperature Tjj ~ 150 MeV |16j . This conjecture was based on his Sta- 
tistical Bootstrap Model. Note, that it was in contradiction to the Fermi's 
model in which the temperature of hadronic matter created at the early 
stage of collisions increases monotonically with collision energy and it is 
unlimited. 

Figure 0] sketches a brief history of pioneering ideas and models on sta- 
tistical hadron production. 

The statistical and hydro dynamical models predict an approximately 
exponential form of particle transverse mass spectra, provided collective 
flow of matter developed in the course of expansion is small. 



Pioneering ideas/models: f(m ) — g" m j/T 

-1950: E. Fermi 

statistical hadron production: T = T~ 

-1951: 1. Pomeranchuk 

freeze-out at T =T FO •« m i 

-1953: L. □. Landau 

hydrodynamical expansion from T to T 
T = f(m, v , T ) 

T FO 

-1965: R. Hagedorn 

statistical hadron production at T = T » ISO MeV 




Fig. 4. Pioneering models and ideas on statistical hadron production. 
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6. S-matrix theory 

The scattering (S-)matrix theory, which relates the initial state and the 
final state of a scattering process, was initiated by Wheeler and Heisen- 
berg in the 30-40s. It led to the development of the Regge theory studying 
analytic properties of the scattering amplitude treated as a function of an- 
gular momentum, and it finally resulted in the string theory pioneered by 
Veneziano. These theories, formulated without specifying elementary par- 
ticles, do not refer to a space-time structure of interaction processes. They 
assume that all particles are bound states lying on Regge trajectories and 
scatter self-consistently. 

Within the spirit of these theories the Wounded Nucleon Model (WNM) 
was proposed by Bialas, Bleszynski and Czyz [10] in 1976. It assumes that 
particle production in nucleon-nucleus and nucleus-nucleus collisions is an 
incoherent superposition of particle production from wounded nucleons, i.e. 
nucleons which interacted inelastically and whose number is calculated using 
the Glauber approach. Up to now, predictions of the WNM model have 
remained an important baseline for interpretation of experimental data. In 
particular, in the case of mean hadron multiplicities the WNM predicts: 

(n AB ) = {w A b)/2 ■ (tinn), (1) 

where (uab) an d ( n NN) ar e mean hadron multiplicities in A+B collisions 
and nucleon-nucleon interactions, respectively, whereas wab is a mean num- 
ber of wounded nucleons in A+B collisions. Note, that the above WNM 
prediction resembles the corresponding prediction of statistical models in 
the grand canonical approximation (thermodynamical models) providing 
the number of wounded nucleons is replaced by a system volume. This ex- 
plains an approximate validity of the WNM for the yield of pions, the most 
popular hadrons. However, the model significantly fails already for yields 
of strange hadrons. Figure [5] presents milestones of the S-matrix era. 

Pioneering ideas/models: 

-1941: W. Heisenberg 

S-matrix theory as a theory of particle interactions 

= 1960: T. Regge + G. Chew, S. Frautschi, J. Collins 
Regge theory 

= 1970: G. Veneziano, S. Mandelstam 
string model 

-1976: A. Bialas, M. Bleszynski, W. Czyz < fg > = <w >/2o<|SJ > 
wounded nucleon model AB AB NN 



Fig. 5. Milestones of the S-matrix era. 
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7. Discoveries of quarks and gluons 

The quark model of hadron classification proposed by Gell-Mann and 
Zweig in 1964 starts a 15 years-long term in which sub-hadronic particles, 
quarks and gluons, were discovered and a theory of their interactions, quan- 
tum chromo dynamics was established. A brief history of this term is shown 
in Fig. [6) In parallel, conjectures were formulated concerning the existence 
and properties of matter consisting of sub-hadronic particles [TTJ [12] , soon 
called a quark-gluon plasma and studied in detail within the QCD [13]. 
The first QCD-inspired estimate of the transition temperature to QGP 
gave T c ~ 500 MeV [H] 

Figure [6] presents a brief history of discoveries of quarks and gluons. 

Many physicists started to speculate that the QGP can be formed in 
nucleus-nucleus collisions at high energies and thus it may be discovered in 
laboratory experiments. 



Pioneering ideas/experiments: 



icatu 




-1964: M. Gell-Mann, G. Zweig 

quark model of hadron classification 

-1965: D. Ivanenko, D. Kurdgelaidze 

quark matter in superdense star cores 

-1968: SLAC experiments: deep inelastic scattering 
discovery of pat tons (now q, q and g) 

-1972: M. Gell-Mann, H. Fritzsch, D. Gross, F. Wilczek, D. Politzer 

qunatum chromodynamics as theory of strong interactions 

= 1975: E. Shuryak 

QCD quark-gluon plasma (T c » 500 MeV) 

-1979: experiments at DESY: three-jet events 
discovery of gluons 



Fig. 6. A brief history of discoveries of quarks and gluons. 
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8. Statistical QGP hadronization and statistical parton 

production 

Questions concerning QGP properties and properties of its transition to 
matter consisting of hadrons were considered since the late 70s. Cabibbo, 
Parisi \15)/ . Hagedorn and Rafelski }16\/ suggested that the upper limit of the 
hadron temperature, the Hagedorn temperature, is the transition temperature 
to the QGP. Furthermore, Rafelski Ji7| / and the collaborators introduced the 
statistical approach to the QGP hadronization. It predicted that the resulting 
system of hadrons is in incomplete equilibrium. The deviations from the 
equilibrium can be traced back to the QGP properties Jiff) /. 

In the mid-90s the Statistical Model of the Early Stage (SMES) was for- 
mulated [19] as an extension of Fermi's statistical model of hadron produc- 
tion. It assumes a statistical production of confined matter at low collision 
energies (energy densities) and a statistical QGP creation at high collision 
energies (energy densities) . The model predicts a rapid change of the colli- 
sion energy dependence of hadron production properties, that are sensitive 
to QGP, as a signal of a transition to quark-gluon plasma (the onset of de- 
confinement) in nucleus-nucleus collisions. The onset energy was estimated 
to be located in the CERN SPS energy range. 

Clearly, the QGP hypothesis and the SMES model removed the contra- 
diction between Fermi's and Hagedorn's statistical approaches (see Sect. 5). 
Namely, the early stage temperature of strongly interacting matter is un- 
limited and increases monotonically with collisions energy, whereas there is 
a maximum temperature of hadron gas, Tq = Tjj ~ 150 MeV, above which 
strongly interacting matter is in a quark-gluon plasma phase. Figure [7] 
summarizes the developments discussed above. 

Pioneering ideas/models: 

~1975: N. Cabibbo, Q. Parisi, 
1978: R. Hagedorn, J. Rafelski 
T c = T H w 150 MeV 

-1991: J. Rafelski 

statistical QGP hadronization 

*1995: M.G.. M. Gorenstein 

statistical production of 
partons at T > T c and of 
hadrons at T < T c 

Fig. 7. Statistical QGP hadronization and statistical parton production. In photo 
from left to right are Johann Rafelski, Tatiana Faberge (former Cern theory group 
Secretary and owner of the rights to Faberge eggs) and Rolf Hagedorn, Divonne- 
les-Bains, 1994. 
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9. QCD based and QCD inspired models 

Attempts to derive from the QCD precise quantitative predictions for 
multi-particle production in high energy collisions have not been success- 
ful [20]. Predictions of QCD-based and QCD-inspired models suffer from 
large quantitative and large qualitative uncertainties, respectively. 

The power-law dependence of transverse momentum spectra, da/dpr ~ 
at high pt (in the perturbative domain) was predicted by Field and 
Feynman in 1977 based on the asymptotic freedom of the QCD. In fact, 
data at px > 2 — 3 GeV follow the power-law dependence, but with the 
power —8 [21] instead of —4. Parton cascade and hadronization models 
were developed [22], however their validity is in question due to statistical 
and hydrodynamical features of experimental data[9], which are difficult to 
explain within the dynamical, QCD-inspired approaches [23] . 

The most famous QCD-inspired models of QGP signals in nucleus- 
nucleus collisions, strangeness enhancement \2J$ and J /ip suppression \25^ . 
were proposed in 1980s by Rafelski, Miiller and Matsui, Satz, respectively. 
They motivated precision systematic measurements. Nevertheless, they did 
not lead to definite conclusions on the QGP creation, because alternative, 
QCD-inspired and/or statistical explanations of data were more successful. 
Figure [8] summarizes the developments discussed above. 



Pioneering ideas/models: f(|TI ) ~ IT1 P 

-1977: R. Field, R. Feynman 

pQCD-based model of high p T phenomena 

~1980: J. Rafelski, B. Mueller, T. Matsui, H. Satz 
QCD-inspired models of QGP signals, 
strangeness enhancement and J/iy suppression 

-1991: K. Geiger, B. Mueller, J. Ellis 

QCD-inspired parton cascade and hadronization model 



Fig. 8. Pioneering QCD-based and QCD-inspired ideas and models. 
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10. Discoveries of strongly interacting matter and its phase 

transition 

Systematic data from experiments at the CERN SPS and LHC and at 
the BNL AGS and RHIC clearly indicate that a system of strongly inter- 
acting particles created in heavy collisions at high energies is close to, at 
least local, equilibrium. At freeze-out the system occupies a volume which 
is much larger than a volume of individual hadron. The latter conclusion 
is based on the failure of the WNM and the success of statistical |26] and 
hydrodynamical models [9]. Thus, one concludes that strongly interacting 
matter is created in heavy ion collisions. 

The phase transition of strongly interacting matter to the QGP was 
discovered within the energy scan program of NA49 at the CERN SPS |27( 
128] . The program was motivated by the predictions of the SMES model. The 
discovery is based on the observation that several basic hadron production 
properties measured in heavy ion collisions rapidly change their dependence 
on collisions energy in a common energy domain |29j . 

Figure [8] summarizes the developments discussed above. 

Pioneering ideas/experiments: 

-1980/00: AGS/SPS/RHIC experiments with heavy ions 
discovery of strongly interacting matter 
(large volume, in ^equilibrium) 



-2000; NA49 at the CERN SPS 

discovery of phase transition of 
strongly interacting matter 
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Fig. 9. Discoveries of strongly interacting matter and its phase transition. The horn 
(left) and step (right) structures in energy dependence of the K + /ir + ratio and the 
inverse slope parameter of K~ mx spectra signal the onset of deconfinement located 
at the low CERN SPS energies. 
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